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Abstract 

The effect of the Hah force on the pinning of vortices in type II superconduc- 
tors is considered. A field theoretic formulation of the pinning problem allows 
a non-perturbative treatment of the influence of quenched disorder. A self- 
consistent theory is constructed using the diagrammatic functional method 
for the effective action, and an expression for the pinning force for indepen- 
dent vortices as well as vortex lattices is obtained. We find that the pinning 
force for a single vortex is suppressed by the Hall force at low temperatures 
while it is increased at high temperatures. The effect of the Hall force is 
more pronounced on a single vortex than on a vortex lattice. The results of 
the self-consistent theory are shown to be in good agreement with numerical 
simulations. 

PACS numbers: 74.60.Ge, 05.40+j, 03.65.Db 
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Typeset using REVTgX 



The advent of high temperature superconductors has led to a renewed interest in vortex 
dynamics. We shall consider the influence of quenched disorder on the vortex dynamics 
in type II superconductors in the presence of a Hall force. The description of the vortex 
dynamics will be based on the phenomenological Langevin equationSS 

R' 

= auRtX n- Vy(R + URt)+FRt + eRt, (1) 

where UR,t is the displacement at time t of the vortex which initially has equilibrium position 
R, rj is the friction coefficient, and m is a possible mass (per unit length) of the vortex. The 
dynamic matrix, $rr', of the hexagonal Abrikosov vortex lattice describes the interaction 
between the vortices in the harmonic approximation. Having a thin superconducting film 
in mind the system is two-dimensional (normal to fi) and the dynamic matrix is specified 
within the continuum theory of elastic mediai by the compression modulus, cn, and the 
shear modulus, Cqq, 

^ ^ 00 / cnql + CrnQy (cn - Cee)qxqy \ 
B \ (cii - C6e)qxqy Ceeq^ + Cug^ j ' 

where (po/B is equal to the area of the unit cell of the vortex lattice, and (pQ = h/2e is the fiux 
quantum. The force (per unit length) on the right hand side of eq. (|Tp consists of the Hall 
force characterized by the parameter a, and Fr^ = (j)o}{Tl, t) x n is the Lorentz force due to 
the transport current density j, and the thermal white noise stochastic force, ^r^, is specified 
according to the fiuctuation-dissipation theorem (^r(^r/(/) = 2r]kBT6(t — t')'^a/3'^RR', and V 
is the pinning potential due to quenched disorder. The pinning is described by a Gaussian 
distributed stochastic potential with zero mean, and thus characterized by its correlation 
function (where now the brackets denote averaging with respect to the quenched disorder) 
{V{x.)V{x.')) = z/(x — x') = z/o/(47ra^) exp(— |x — x'p/(4a^)), taken to be a Gaussian function 
with range a and strength uq. 

Upon averaging with respect to the quenched disorder the average restoring force, F/? = 
^ Er' '^RR'((uR't)), of the lattice vanishes. On the average, corresponding to the lattice 
reaching a steady state velocity v = ((li)), there will be a balance, F + Fj + + Fp = 0, 
between the Lorentz force, F, the friction force, Ff = — r^v, the Hall force, Fh = av x n, 
and the pinning force, Fp = —{{W)). The pinning force is due to time-reversal symmetry 
invariant under reversal of the direction of the magnetic field, and is therefore antiparallel to 
the velocity.0 Thus, the pinning yields a renormalization of the friction coefficient in terms 
of a velocity dependent effective friction coefficient, Ff + Fp = —ries{v)\, which reduces 
in the absence of disorder to the bare friction coefficient t], and has previously only been 
determined to lowest order in the disorder.! The relationship between the average vortex 
velocity and the induced electric field, E = v x B, leads to the expressions for the resistivity 
tensor and Hall angle 

p = -J^ ^ ( '^off \ Q ^ arctan — . (3) 

vis + Ves J rj^s 

The average vortex motion is conveniently described by reformulating the stochastic 
Langevin problem in terms of a path integral. The probability functional for a realization 
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{uRt}R of the motion of the vortex lattice may be expressed, using the equation of motion, 
through a functional integral over a set of auxiliary variables {uRt}R, and we are led to 
consider the generating functionali'i 

Z[F, J] = /n ^5uR, / n Vun^, Je^^I"'"], (4) 

R R' 

where in the action, S[u, u] = u{{D^)^^u + F — VV^ + + Ju> have introduced a source 
field J coupling to the vortex positions u, and used matrix notation in order to suppress 
the integrations over time and summations over vortex positions and Cartesian indices. The 
retarded Green's operator is given by —{D^)~^u = muju + V^m + J^-R' '^'RR'UR't + anx ur^, 
and its Fourier transform is 

(D^)-^ = ( "^^^ + ^^^ ^ _ $ (^) 

In order to immediately be able to perform the average with respect to both the Langevin 
noise and the disorder, we have chosen a non-zero mass, m 7^ 0, leaving the Jacobian, J', 
an irrelevant constantilll (in final expressions the mass can be set to zero, and will in fact 
for the values chosen not affect the obtained numerical results) and we obtain the averaged 
functional 

Z[f] = {{Z))= [V(j)e'^^'^^+'f^ (6) 



which generates, for example, the average position and correlations 



^((URt)) 



6J 



Rt 



, ((URtUR/t/)) 



(7) 
j=o 



We have introduced the notation (f) = (u, u) and / = (F, J), and the action, S = Sq + Sy, 
consists of a quadratic term, So[(f)] = (f)D~^(f)/2, where the matrix in addition is a matrix 
in Cartesian indices, and time and vortex positions (5rr/ = 5rr' S{t — t')) 

and a term originating from the disorder 

iSv[(f>] = ^ E Jdtjdt' u^tdadpu{ujit - uji>t')un't'- 



RR'a/3 ' 



(9) 



This reformulation of the stochastic problem in terms of a field theory is equivalent to the 
formalism of Martin, Siggia and Rose,! as noted previously.! 

Our aim is to express the effective action in terms of all two-particle irreducible vacuum 
diagrams, and we therefore add a two-particle source term to the generating functional 

= /p0e*^[*l+*^'^+^*^<^. (10) 
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The generator of connected Green's functions, i^] = —i\nZ[f,K], has accordingly 

derivatives given by (the bar consequently denotes the average with respect to the action 
S[0] + /0 + 0ir0/2) 



6W 



SW 



5/, 



a 



"-"■RtR't' 



\i4>m <Pwt' + iGa^iKt, R't')), 



(11) 



where G is the full connected Green's function of the theory. The quantity of interest is the 
effective action r[0, G] = W[f, K] — fcj) — <pK(p/2 — iGK/2^ the Legendre transform which 
satisfies the equations 



5V_ 

^5 



Jg 



-K. 



(12) 



In the physical problem of interest the sources K and J are absent, K = {] and J = 0, and 
the full matrix Green's function has, due to the normalization of the generating functional, 
Z[F, J = 0, is: = 0] = 1, the structure 



G, 




' (((5m" bu^))\ 
^{{Su^ Su^)) {{Su" Su^ j) ) 



(13) 



where 6u = u— ((u)) and 6u = u— ((u)). The retarded Green's function gives the linear 
response to the force F/3, and is the correlation function (both matrices in Cartesian 
indices as indicated). 

According to Cornwall et a/.,0 the effective action can be written on the form r[(/), G] = 
S[^] + ^Tt{{Ds^^ lnD-i)G - 1) - i ln(e*^'"O^P\ where Dg^ = 6^S[^]/5^S^, and ^int[^,0] 
is the part of S[4> + ip] which is higher than second order in ip in the expansion around 0, 
and Tr denotes the trace over all variables. The superscript "2PI" on the last term indicates 
that only the two-particle irreducible vacuum diagrams should be included in the interaction 
part of the effective action, and the subscript that propagator lines represent G, i.e., the 
brackets with subscript G denote the average {F[iP])g = (detG)-i/2/p^ e'^^~'^/^F[iP], for 
an arbitrary functional F. We now expand the exponential and keep only the first order 
term in S'int and obtain 



(14) 



For the physical problem of interest the two particle source, K, vanishes, and the last of 
the equations in ([T2|) therefore yields the Dyson equation, G^^ = D^^ — G], with the 
matrix self-energy given by 




2i 



.6{Sy[^ + ^P])f' 



SGij 



(15) 



K=0,J=O 



The Dyson equation and eq.(^) constitute a set of self-consistent equations for the 
Green's functions and the self-energies. The average field occurring in eq.(|T5[) is given 
hj (j) = (((u)), ((uRt))) = (0,vt), as the expectation value of the auxiliary field van- 
ishes, ((u)) = —i Z^^6Z/6F\j^Q j^^Q = 0, due to the normalization of the generating 
functional. The matrix self-energy has two independent components, and (as 
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S^„(Rt,Rr) = [Sf^(R't',Rt)]*), and for vortices we have according to eq. (|T5D 
= -<=o,a.=o> where af^(RtR't') = 1/iVEk ^^(k)A;,A;^(kG«(RtR't')k)e^'^^ and 
Sf^(Rt, H't') = —i/NJ2k v{)i)ka_kpe^'P^ . The influence of thermal and disorder induced fluc- 
tuations is described by the phase, ip\^ = ikMk + k ■ (R — R' + ur^ — ur/^/), specified by 
the Cartesian matrix Map(Rt,K't') = z(Gf^(Rt,Rt) - ^^^(Rt, R't')). Using the Langevin 
equation and the first equation in (|12D we obtain for the pinning force 

F =- 



J2 Jdt' J2ku{k){kG^{mR't')k)e'^^. (16) 



We first consider the case of non-interacting vortices. This is appropriate for low magnetic 
fields where the vortices are so widely separated that the interaction between them can be 
neglected. We have solved the above set of self-consistent equations by numerical iteration. 
In fig. Il| the resulting pinning force as a function of the velocity is shown for a set of different 
strengths of the Hall force in the low temperature regime, i.e., T ^ z/q /{ksa). The Hall 
force is seen to reduce the pinning force in this temperature regime except, of course, at low 
velocities. The high velocity behavior, v ^ \f^l{j]Q^^-, can be compared with the second 
order perturbation expression, which is obtained by replacing the full Green's function in 
eq. (0) by the free Green's function, and omitting M in the exponent (the mass term can 
be neglected assuming m <^ rfa? /^/U^) 

According to fig. ^ there is good agreement between the self-consistent and perturbation 
theory in the reduction of the pinning force due to the Hall force at high velocities. 

At high temperatures, T u^^"^ /{ksa), and moderate velocities, v < y^/{ria'^), the 
Hall force has the opposite effect on the pinning force. According to eq. (|l^) we obtain (for 

- -SnviksTya- ^^^^ 

In this high-temperature limit (which can be realized in high temperature superconductors) 
we observe that the self-consistent theory yields a pinning force that has a linear velocity 
dependence and that the Hall force yields an increase of the pinning force, as shown in the 
inset in fig. 0. 

In order to fully test the validity of the self-consistent theory its results are also compared 
to numerical simulations as shown in fig. ^ The agreement between the self-consistent 
theory and the simulations is good except around the maximum of the pinning force. In this 
region the relative velocity fiuctuations are large and the self-consistent theory predicts that 
the relative fiuctuations are diverging at zero velocity even at T = 0. The self-consistent 
equations (as well as the numerical simulations) can therefore be expected to yield the largest 
errors at low velocities. 

The Hall angle is from the self-consistent theory found to increase monotonically from 
zero at low velocities to the disorder independent value 9o = arctan(a/r7) at high velocities, 
as shown in fig. El for the single vortex case. The agreement between the self-consistent 
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theory and the numerical simulations is seen to be good, testifying to the validity of the 
approximation made in eq. (14). As shown in fig. ^we find that increasing the temperature 
increases the Hall angle at low velocities and that this feature vanishes at high velocities. 

Finally we consider a vortex lattice treating the interaction between the vortices in the 
harmonic approximation. The pinning force obtained from the self-consistent theory for the 
case of zero temperature is shown in fig. |. As expected there is no influence of the Hall force 
on the pinning force at low velocities, but we find a suppression at intermediate velocities, 
and at high velocities, v ^ Cua/rj, we recover the high velocity limit of the single vortex 
result, i.e., eq. (|T7p. By comparison of fig. |I| and fig. ^ we find that the Hall force has 
a much weaker influence at intermediate velocities on the pinning of an interacting vortex 
lattice than on a system of non-interacting vortices. The influence of the Hall force on the 
pinning force is more pronounced for a stiff than a soft lattice as seen from the inset in fig. 
^, and is similarly reflected in the Hall angle dependence on the stiffness of the lattice as 
seen from the inset in fig. ^ the stiffest lattice has the greatest Hall angle. 

A possible experimental verification of the obtained results would be to measure the Hall 
angle and pinning force of a type II superconductor, and thereby obtain the value of a of the 
particular material according to eq. (^. The parameters characterizing the disorder, a and 
1^0, may, e.g., be determined by both measuring the velocity dependence of the pinning force 
at high vortex velocities and at high temperature at moderate velocities. The self-consistent 
theory can then be compared to the experimental results for pinning forces and Hall angles 
using the experimentally obtained parameters as input. 

In conclusion, we have studied analytically as well as through simulations the vortex 
dynamics in type II superconductors in the presence of a Hall force and quenched disorder. 
For the case of a single vortex we find that the Hall force reduces the pinning force in the 
high-velocity regime where the influence of fluctuations is negligible and the only effect of 
the Hall force is through the response function. The situation at high temperatures is the 
opposite since then the thermal fluctuations are dominating over the influence through the 
response function, and the Hall force thus increases the pinning force because it suppresses 
the fluctuations. The influence of the Hall force on a vortex lattice is found to be weaker 
than on a single vortex. 
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FIGURES 
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FIG. 1. Pinning force (in units of fg^^a^^) as a function of velocity (in units of ?7~^a^^fg^^) 
for a single vortex for various strengths of the Hall force. The curves correspond to 
a/rj = 0,0.2,0.4,0.6,0.8,1, where the uppermost curve corresponds to a = 0. The mass is 
m = ^^"^ and the temperature is T = 0.1z^y^/(/csa). Inset: Pinning force (in units 

of 10~^i'g^^a~^) as a function of velocity (in units of ?7~^a~^fg'^^) according to the self-consistent 
theory at high temperature, fesTa/z/y^ = 10. The upper curve corresponds to a = ??, the lower to 
a = 0. The mass is m = O.Olrj^a^UQ 
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FIG. 2. Comparison of the simulation results for the pinning force and the results of the 
self-consistent and second order perturbation theory for a single vortex for the case of no Hall force 
(a = 0) and a moderately strong Hall force {a = rj). The solid line represents the self-consistent 
result and the crosses the simulation result while the uppermost dashed-dotted line represents 
the perturbation theory result, all for the case a = 0. The dashed line and the plus symbols 
represent the self-consistent and simulation results, while the lowest dashed-dotted line represents 
the perturbation theory, all for a = rj. The mass is m = O.It^^o^z/q ^''^ and the temperature is 
T = Q.lvy^ /{kBo). The units of the pinning force and velocity are chosen as in fig. |^. 
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FIG. 3. Hall angle for a single vortex as a function of velocity. The curves represent the 

— 1/2 

self-consistent results for the temperatures ksTav^ = 0,0.1,1, where the uppermost curve 

corresponds to the highest temperature. The plus symbols represent the simulation result for 

ksTauQ = 0.1. The parameter a/r/ is unity and the mass is 0.177^0^2^^^. Inset: Hall angle for 

1 12 

a vortex lattice as a function of velocity in descending order of lattice stiffnesses Ac^q = 200z/q , 
lOOfQ^^, SOt'g^^. The unit of the velocity occurring in the figures is chosen as in fig. |l[ 




FIG. 4. Pinning force (in units of Uq A^'^/'^) as a function of velocity (in units of tj^^i'q A^^^^) 

for a vortex lattice of size 16 x 16. The range of the disorder correlator, a, is chosen to be 0.1^^^^, 

where A is the unit cell area. The solid curve corresponds to a = while the dashed corresponds to 

a = rj. The temperature and mass are both set to zero. The elastic constants are Acn = lO^fg''^ 
1/2 

and Acqq = IOO^q . Inset: Pinning force as a function of velocity for a = and a = rj, respectively. 

1/2 

Here Acqq = SOOt'g and the other parameters are chosen as above. 
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